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Abstract 

We study spectral features in the gamma-ray emission from dark matter (DM) annihi¬ 
lation in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), with either 
neutralino or right-handed (RH) sneutrino DM. We perform a series of scans over the 
NMSSM parameter space, compute the DM annihilation cross section into two photons 
and the contribution of box-shaped features, and compare them with the limits derived 
from the Fermi-LAT search for gamma-ray lines using the latest Pass 8 data. We imple¬ 
ment the LHC bounds on the Higgs sector and on the masses of supersymmetric particles 
as well as the constraints on low-energy observables. We also consider the recent upper 
limits from the Fermi-LAT satellite on the continuum gamma-ray emission from dwarf 
spheroidal galaxies (dSphs). We show that in the case of the RH sneutrino the constraint 
on gamma-ray spectral features can be more stringent than the dSph bounds. This is 
due to the Breit-Wigner enhancement near the ubiquitous resonances with a CP even 
Higgs and the contribution of scalar and pseudoscalar Higgs final states to box-shaped 
features. By contrast, for neutralino DM, the di-photon final state is only enhanced in 
the resonance with a Z boson and box-shaped features are even more suppressed. There¬ 
fore, the observation of spectral features could constitute a discriminating factor between 
both models. In addition, we compare our results with direct DM searches, including the 
SuperCDMS and LUX limits on the elastic DM-nucleus scattering cross section and show 
that some of these scenarios would be accessible to next generation experiments. Thus, 
our findings strengthen the idea of complementarity among distinct DM search strategies. 
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1 Introduction 


The existence of a vast amount of dark matter (DM) in the Universe is supported by many 
different and independent observations. The latest measurements suggest that approximately 
a 27% of the Universe energy density is in form of a new type of non-baryonic cold DM [T]. 
Given that the Standard Model (SM) of particle physics does not contain any viable candidate 
to account for it, DM can be regarded as one of the clearest hints of new physics. Among 
the many diverse proposals within particle physics, Weakly Interacting Massive Particles 
(WIMPs), provide a very well-motivated paradigm, since these particles can be thermally 
produced in the early Universe in a sufficient amount to reproduce the observed relic abun¬ 
dance. 

The hunt for DM is an extremely active field. These exotic particles could be observed 
through their collisions with nuclei inside direct detection experiments, a challenging tech¬ 
nique that involves the use of very sensitive underground detectors. They could also be 
produced at colliders, such as the LHC, where they are expected to leave an imbalance in the 
transverse energy as they escape the detectors. Finally, the annihilation (or decay) of DM 
particles in the Galactic halo as well as in other astrophysical objects, and the subsequent 
emission of SM particles (gamma-rays, antiparticles or neutrinos) is being searched for by a 
variety of techniques, collectively referred to as indirect detection. The challenge in the latter 
strategy is the discrimination between a potential signal and the astrophysical background, 
which requires a good modelling of all the conventional physical processes that can produce 
cosmic rays and photons. 

Although various potential hints of DM detection have been reported in the last decades, 
both from direct as well as indirect detection experiments, none of these have been conclu¬ 
sively confirmed yet. To their credit, the DAM A/LIBRA annual modulation [2]; the 511 keV 
excess observed by the SPI experiment [3]; the PAMELA and AMS rise in the positron frac¬ 
tion miS]; a potential X-ray line emission near 3.5 keV observed in our Galaxy and certain 
galaxy clusters [HE]; and the apparent excess in Fermi-LAT data from gamma-rays in the 
Galactic Gentre (GG) [SltlTj have triggered an enthusiastic examination of DM models be¬ 
yond the “vanilla” scenario during the past decade. Among these potential signatures, the 
Fermi-LAT low-energy excess in the Galactic Centre Emission (GCE) [TSj stands out, as it 
could be explained in terms of WIMP DM, with a mass in the 6 — 200 GeV range and an 
annihilation cross section compatible with that of thermal relief. In any case, the verification 
of this and other DM hints would need an independent complementary measurement, which 
would allow us to gain further insight into the DM particle nature. 

In this regard, the observation of spectral features in the gamma-ray spectrum would 

^Although an astrophysical explanation in terms of millisecond pulsars, or leptonic cosmic ray outbursts 
seems also plausible m- 
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constitute an excellent discrimination tool, as these features are linked to characteristic prop¬ 
erties of DM models. For example, the energy of monochromatic photons from either the 
annihilation or decay of DM particles could provide a measurement of the DM mass. Simi¬ 
larly, box-shaped features and internal bremsstrahlung processes can be associated to specific 
annihilation channels. Furthermore, since in general, astrophysical processes predict contin¬ 
uous 7 -ray spectra at the GeV energy scale, this observation would be a smoking gun of DM 
detection. 

The Fermi Large Area Telescope (Fermi-LAT) is currently performing an all-sky survey 
in 7 -rays with unprecedented sensitivity in the GeV range. Specially interesting for DM 
searches is the GG region, where a high density of DM particles is expected and therefore 
higher gamma-ray fluxes from their annihilation or decay are predicted, as well as nearby 
dwarf spheroidal galaxies (dSphs), which are DM dominated objects. No DM signal has 
been obtained from the latter, leading to very stringent upper bounds on the current DM 
annihilating cross section in the Galactic halo, {crv)Q, which rules out the “thermal” value 
{{crv)o Ri 3 X 10“^® cm^/s) for WIMP masses below ~ 100 GeV (the precise value depends 
on the main annihilation channel). These constraints are in mild tension with the DM 
interpretation of the GGE that we have briefly mentioned above |20ll21j . 

Regarding the search for features in the gamma-ray spectrum, no significant global excess 
of a spectral line has been found yet. A recent analysis of 5.8 years of data (Pass 8 ) [22] 
has explored the energy range from 200 MeV to 500 GeV for five different regions of interest 
(ROIs), optimised to maximise the signal to noise ratio. A previous study had extended the 
search to energies as low as 100 MeV [23] using the Pass 7 data set. Upper constraints on the 
annihilation cross section of DM particles into monochromatic photons have been derived in 
the range {av)^^ ~ 10“^® — 10“^^ cm^/s, depending on the WIMP mass and the DM density 
profile assumed. These limits begin to explore the expected values for a loop-suppressed 
WIMP annihilation cross section and thus constrain some regions of the parameter space of 
DM models. 

From the perspective of particle models for DM, the annihilation of WIMPs into a pair 
of photons is a loop-suppressed process, hence it is expected to be three or more orders of 
magnitude below the thermal cross section |24H36j . However, there are known mechanisms 
by which the line emission can be enhanced, which were explored in order to account for a 
previous apparent excess at 130 GeV [37IH2] that is no longer statistically significant. For 
instance, resonant annihilation, together with annihilation just above the production thresh¬ 
old of the particles involved in the loop [l3] can significantly increase the rate of this process. 
Furthermore, gamma rays can also be produced in the decay in-flight of certain annihila¬ 
tion final states (such as scalar and pseudoscalar Higgs bosons), giving rise to characteristic 
box-shaped features in the gamma-ray spectrum. Finally, gamma-ray features from internal 
bremsstrahlung |44ll45j . linked to processes in which charged particles mediate DM annihila- 
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tion, can also lead to observable signatures. The theoretical predictions for these processes 
are extremely model dependent and, as such, could be used in order to distinguish between 
different DM scenarios. 

In this article, we work within the context of the Next-to-Minimal Supersymmetric Stan¬ 
dard Model (NMSSM). The NMSSM is a well-motivated extension of the Minimal Supersym¬ 
metric Standard Model in which an extra singlet field, S, is introduced, which mixes with 
the Higgs SU(2) doublet and whose vacuum expectation value after electroweak symmetry 
breaking generates an effective EW-scale /U parameter [46] (see, e.g.. Ref. |47| for a review). 
Among its many virtues, the NMSSM displays a very interesting phenomenology, mainly due 
to its enlarged Higgs sector. The mixing of the Higgs doublet with the new singlet component 
opens the door to very light scalar and pseudoscalar Higgs bosons with interesting prospects 
for collider searches. Moreover, in the NMSSM the mass of the Higgs boson also receives 
new tree-level contributions from the new terms in the superpotential [481149] . which allow 
to increase the Higgs mass and reduce the amount of hne-tuning of the model |50H52j . The 
profound changes in the neutralino and Higgs sectors also have interesting implications for 
neutralino DM and new contributions to both direct and indirect detection. The NMSSM 
can be enlarged with an extra singlet superfield that incorporates right-handed neutrinos 
(and sneutrinos) [531154] in order to accommodate a see-saw mechanism that explains the 
smallness of neutrino masses. The right-handed (RH) sneutrino in the resulting construction 
is a viable DM candidate [55] with interesting phenomenological properties. 

In this paper, we compute the annihilation cross section into two photons and the contri¬ 
bution from box-shaped features for the lightest neutralino and the RH sneutrino in the afore¬ 
mentioned supersymmetric constructions, focusing on the low-mass region below 200 GeV. 
We compare the results with the current upper bounds derived from the search for features in 
the gamma-ray spectrum studied by the Fermi-LAT collaboration. We investigate the effect 
of these constraints and compare them with Fermi-LAT limits on the continuum emission of 
gamma-rays from dSphs, as well as direct detection bounds on the elastic scattering cross 
section. It is important to emphasize that we do not make any attempt to fit the GCE in 
this paper. We observe that for certain regions of the parameter space with RH sneutrino 
DM, the search for these spectral features can be more constraining than the existing bounds 
from dSphs, and even than direct detection limits. This rarely happens for the case of the 
neutralino, unless the narrow resonant annihilation condition with the Z boson is fulfilled. 
These results emphasize the relevance of this technique in order to discriminate among DM 
candidates. 

This article is organised as follows. In SectionOwe investigate the NMSSM with neutralino 
DM and in Section[3]we address the case of the RH sneutrino in an extended NMSSM. Finally, 
the conclusions are presented in Section [H 
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2 Spectral features from neutralino annihilation in the NMSSM 


The modifications in the neutralino sector due to the inclusion of a singlino component and 
the presence of new annihilation channels have profound consequences for DM searches, which 
have been extensively discussed in the literature [56H6nj . For example, it has been shown 
that the NMSSM can accommodate low-mass neutralino DM |31ll61f[67] . Several studies have 
recently pointed out that these neutralinos can also explain the GCE [681172] when they are 
either singlino-Higgsino or bino-Higgsino admixtures [ 68 ]. 

The monochromatic gamma-ray spectrum from neutralino annihilation in the context of 
the NMSSM has been previously studied in Refs. [31113511731175] . In general, one expects that 
the cross section of the process XiXi 77 is loop suppressed with respect to the thermal 
annihilation cross section. In Ref. m, the authors showed that for neutralinos lighter than 
100 GeV in the NMSSM, the cross section for the 77 final state is typically of the order of 
10“^^ cm^ s“^ for either bino-like, singlino-like or mixed neutralinos. This value can increase 
several orders of magnitude when neutralinos have a resonant annihilation mediated by a 
very light CP-odd Higgs with a fermion loop, reaching {crv)^^ ks 10 “^^ cm^s“^ [3TH731175] . 

In this section, we update these results, incorporating the LHC constraints on the Higgs 
sector. They set stringent bounds on the Higgs couplings and the invisible branching ratio 
and therefore also limit the availability of very light scalars and pseudoscalars, needed in 
order for the light neutralinos to fulfil the relic density constraint [6311641176] . In addition, 
we take into account the most recent bounds from direct and indirect DM searches. To this 
aim, we first compute the monochromatic gamma-ray spectrum for a series of scans over the 
parameter space and determine the observability of this signal. Then, we compare the results 
with other direct and indirect search strategies. 

When computing the continuum contribution to the gamma ray spectrum, we must also 
take into account virtual internal bremsstrahlung (VIB) processes that can take place in t- 
channel exchange of a charged particle. In the case of the neutralino, this can correspond to 
annihilation channels into lepton pairs (mediated by sleptons), H^H~ (mediated by charged 
Higgses), or W~^W~ pairs (mediated by charginos). The contribution from VIB can exceed 
that from final state radiation [HITT] and in some cases be more easily observable than 
gamma-ray lines [78] . As it has been shown in Ref. m, the effect of VIB is enhanced when 
the mass of the DM and the mass of the charged mediator are almost degenerate. We will 
later argue that these conditions are not satisfied in our scan. 

2.1 Details of the scan and experimental constraints 

We have carried out a series of scans, designed to explore neutralino masses below 200 GeV. 
The NMSSM input parameters have been dehned at the EW scale according to the ranges 
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in TablelH Fixed values are used for the gluino soft mass, M 3 = 1500 GeV, for the trilinear 
parameters, Ajj = 3700 GeV, Ad = 2000 GeV, and Ae = —1000 GeV, as well as for the 
soft scalar masses of sleptons and squarks, = m^, = 300 GeV and ttIq^ = = 

1500 GeV, respectively, where the index i runs over the three families. The conservative choice 
of the squark masses is motivated by the LHG null results from supersymmetry searches. Also 
note that despite the large trilinear term Ajj, the instability against charge- and/or color¬ 
breaking minima is avoided since the soft mass of the squarks is at the TeV scale [80] . A large 
Ajj helps to increase the loop contributions to the GP-even Higgs mass and thus is useful for 
reproducing the experimental value (this is more relevant in the MSSM than in the NMSSM, 
where extra tree-level contributions are present). Small slepton masses are helpful to enlarge 
the supersymmetric contribution to the muon anomalous magnetic moment (although this 
constraint is not imposed in our analysis). We consider that these scans are sufficient to 
study the phenomenology of monochromatic gamma-ray emission from the annihilation of 
low-mass neutralinoj^. 

The range of variation of our parameters in Tabled) is similar to the one used in other 
scans of low-mass neutralino DM in the NMSSM [631167[ [76]. Some of these analyses have 
used MCMG techniques to explore the NMSSM parameter space. In the case of Ref. [63] 
the authors considered 11 initial parameters. We have focused our scan in the regions where 
more solutions were found in Ref. [63], but adopting fixed values for the scalar masses. Also, 
Ref. [67] extends the range in |Ak| up to 200 GeV, but in Ref. [63] small values of |Ak| are 
preferred in order to have low-mass neutralinos. 

The masses of the supersymmetric particles and low-energy observables have been com¬ 
puted with NMSSMTools 4.1.2 |811183j . whereas the gamma-ray spectrum and the DM relic 
abundance have been computed with micrOMEGAs 3.6.9 [HI]. In order to efficiently explore 
the parameter space, we have linked these codes with MultiNest 3.9 [851187) . which uses a 
likelihood function to generate MCMG and to find regions of the parameter space that max¬ 
imise the likelihood. This function contains uncorrelated Gaussian probability distributions 
for the DM relic abundance, the SM-like Higgs mass, BR(Hs —^ and BR (6 —>■ 57 ), 

centred around the observed values. 

It is important to stress that the scan procedure has not been optimised to extract any 
statistical information on the results and should not be regarded as a complete exploration 
of the vast multidimensional parameter space. In this article, we are only concerned with 
obtaining a representative subset of viable solutions that covers the low-mass range of the 
DM parameter space without addressing the probability of each point to explain the whole 
set of experimental data considered. 

^As we will comment later, we do expect that modifications of some of these input parameters have an 
influence on other aspects of the neutralino phenomenology, such as the contribution of VIB to the continuum 
gamma-ray spectrum. 
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Parameter 

Scan 1 

Scan 2 

Scan 3 

Ml 

[ 1 , 200 ] 

[1,40] 

[ 1 , 200 ] 

M 2 

[ 200 , 1000 ] 

[ 200 , 1000 ] 

[700,1000] 

tan/3 

[4,20] 

[4,20] 

[2,50] 

A 

[ 0 . 1 , 0 . 6 ] 

[ 0 . 1 , 0 . 6 ] 

[ 0 . 001 , 0 . 1 ] 

K 

[ 0 , 0 . 1 ] 

[ 0 , 0 . 1 ] 

[ 0 . 1 , 0 . 6 ] 

Aa 

[500, 5000] 

[500, 5000] 

[500,1100] 

A. 

[-50,50] 

[-30,0] 

[-50,50] 

f^eff 

[110,250] 

[160, 250] 

[200,400] 


Table 1: Input NMSSM parameters for the series of scans used in the neutralino case. Masses and 
trilinear parameters are given in GeV. All parameters are defined at the EW scale. 


For each point of the parameter space, we have computed the full supersymmetric spec¬ 
trum and imposed the most recent experimental constraints from collider searches. For the 
sleptons and charginos, we include the lower limits on their masses from LEP through the 
micrOMEGAs 3.6.9 code. For all the squarks, we apply a lower limit on the mass of 1.5 TeV 
independent on the gluino mass, which is in agreement with the latest ATLAS results [ 88 j . 
Besides, our choice for the gaugino soft mass parameter, M 3 , ensures that the gluino mass 
is above the current LHC bound. Regarding the Higgs sector, we have ensured that one of 
the CP-even Higgs bosons lies in the mass range 123 GeV< m^o < 128 GeV (that accounts 
for the experimental and theoretical uncertaintjH) and has properties very close to those of 
the SM Higgs. In particular, we have imposed an upper bound on its invisible branching 
fraction, —>■ inv) < 0.27. We also have included the bounds from the rare de¬ 
cays Bs —>■ ■, b ^ s'y and —>■ For further details on how these bounds are 

implemented, see Ref. m- 

We have also implemented the most recent constraints from DM direct detection ex¬ 
periments. In particular, we have made use of the latest results of LUX [92] and Super- 
CDMS [931194] . which are the most constraining experiments in the mass range considered 
along this worl0. We have computed the upper bound for these two experiments in each 
point of the parameter space of the scanned model using the Yellin’s maximum gap method 
at 90% C.L. |96| following the procedure sketched in Ref. m- With the data provided in 
that work for different DM haloes, the spin-independent and spin-dependent contributions 
are tested simultaneously using the expected number of events in the corresponding detec¬ 
tor, thus deriving consistent limits. For consistency with the analysis of gamma-ray lines 

®This range might be enlarged due to two-loop corrections to the Higgs masses in the NMSSM, which can 
be especially important for large values of A |89II90| . 

For smaller masses, the CRESST collaboration provides the most stringent bound on the elastic scattering 
cross section of DM particles up to date | 95 |. 
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of Ref. [22], the direct detection bounds have been extracted assuming 
density, 0.4 GeV/cm^, and the same choices of DM halo density profile: 


the same local DM 


We have set an upper bound on the DM relic abundance, < 0.13, consistent with 
the latest Planck results [T]. We have also considered the possibility that supersymmetric 
DM only contribute to a fraction of the total relic density, and set for concreteness a lower 
bound on the relic abundance, 0.001 < D/i^. The fractional density, ^ = min[l, Dfi^/O.ll], 
has been introduced to account for the reduction in the rates for direct and indirect searches 
(assuming that the DM candidate is present in the DM halo in the same proportion as in the 
Universe n. 


Finally, regarding indirect detection, we have also incorporated the Fermi-LAT bounds 
on dSphs, we remind the reader that we do not attempt to fit the GCE. The Fermi-LAT 
collaboration has performed an analysis of the gamma-ray emission from 25 dSphs using six 
years of data |2T|. The absence of a signal can be interpreted as constraints on the annihilation 
cross section of DM particles. It is customary to assume annihilation into pure SM channels 
in the calculation of these bounds. Nevertheless, the occurrence of non standard annihilation 
final states in our DM scenarios prevents us from using these results directly. Instead, we have 
extracted independent 95% GL upper bounds on C^(cru)o for the six more constraining dwarf 
galaxies, i.e. the ones with the largest J-factors [99| (Goma Berenices, Draco, Segue I, Ursa 
Major II, Ursa Minor and Willman I). We have derived the upper limits from the gamma-ray 
flux predicted by the specific DM model for each dSph, using the delta-log-likelihood method, 
the mean values of the J-factors and the bin-by-bin likelihood functions provided in Ref. [99j . 
without taking into account the uncertainties on the kinematically measured J-factors. Then, 
we have applied the most restrictive of these limits to our data. We have checked that this 
procedure leads to slightly less stringent bounds than the combined limit from the Fermi-LAT 
collaboration (by a factor smaller than 1.5 in the whole mass range), when applied to the 
region of DM masses from 10 to 100 GeV with pure annihilation channels. To incorporate 
the latest results from dSphs |21] (Pass 8 analysis), we have estimated the improvement of 
this data set with respect to the previous analysis [99| in a factor 4 on for all DM 


masses. 




Figure 1: Lightest pseudoscalar mass versus the lightest neutralino mass for all the points that 
satisfy the experimental constraints. The different colours denote the main annihilation channel in 
the Galactic halo, as indicated in the legend. 

2.2 Results 

Figured] shows the mass of the lightest pseudoscalar as a function of the mass of the lightest 
neutralino for all the points in the scan that satisfy the different experimental constraints. 
Not surprisingly, the upper bound on the neutralino relic density is responsible for removing 
extensive regions of the parameter space. As we can observe from this figure, the vast 
majority of the solutions found satisfy the condition for resonant annihilation with either the 
lightest pseudoscalar (when m^o ~ 1/2m^o, as indicated by the dashed line), or with the Z 
boson, or with the SM-like CP-even Higgs (the two narrow vertical bands for m^o ~ 1/2 
and m^o ~ respectively). The correct relic density can also be obtained when 

the annihilation into a pair of very light pseudoscalars is kinematically allowed (below the 
solid line that corresponds to m~o = rn^o). Some points in our scan appear to cluster 
around the line with m^o = m^o, but these in fact correspond to cases in which the resonant 
condition with the second-lightest Higgs is satisfied, when is heavier than the SM Higgs. 
These conditions are extremely difficult to achieve and have been widely discussed in the 
literature [SI1E31IM1EZ1IZS]. Given the low-mass range considered in this analysis, and our 

®Note that the profiles (both Einasto and NEW) considered in Ref. are slightly different from those of 
Ref. |98| . used to extract the speed distributions in Ref. m- However, we have checked that these differences 
only lead to small deviations in the expected number of events. 

®This assumption is valid if another cold DM candidate, e.g., the axion, accounts for the missing relic 
density. Changes in the thermal history of the Universe or the presence of DM candidates with different 
interaction properties, such as self-interacting DM, might modify this hypothesis. 
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choice of the soft masses for sleptons and squarks, no coannihilation effects are present (only 
a few points present coannihilation with the second-lightest neutralino). 

In Figure [21 we represent the thermally averaged neutralino annihilation cross section into 
two photons in the Galactic halo and compare it with the Fermi-LAT upper bounds, derived 
for the Einasto and NFWc DM density profiles in their corresponding optimised ROIs [22j . 
Although the results from Ref. [23] extend the bounds to lower masses for the Einasto profile, 
they do not exclude any new points in our scan. The Eermi-LAT bounds are shown as solid 
black lines with their corresponding 68% and 95%CL expectations (black dotted and thin 
black lines respectively). Note that our theoretical predictions are in general very far from 
the current constraints. Nevertheless, large values of can be obtained in the vicinity 

of the resonances due to the well-known Breit-Wigner enhancement, which has been studied 
in various models |421I651[T00[ - I102j . The most obvious resonance is the one with the Z boson, 
which appears as a prominent peak for m^o ~ 45 GeV. The resonance with the lightest 
pseudoscalar or with a GP-even Higgs is not restricted to a single value of neutralino mass; 
due to the variation in the Higgs mass we can observe points with an increased 
for the whole mass range (the peak for m^o ~ 62 GeV corresponds to the resonance with 
the SM-like Higgs). Notice that the conditions for a Breit-Wigner enhancement are more 
easily satished in the Z boson resonance, given its larger decay width, whereas for Higgs 
boson resonances this only happens for an extremely small range of masses. In the neutralino 
case we must therefore conclude that the observation of gamma-ray lines by the Eermi-LAT 
satellite is very unlikely, and subject to exceptional conditions. 

A small fraction of our solutions, with rr final states, feature stau exchange, and therefore 
potential contribution from VIB. This occurs for some of the red points in Eigure 2 with 
~ 45 GeV and ~ 60 GeV. In Ref. [37], it is shown that when = (ruf/m^)^ is 
large, the VIB contribution scales as In our scan, the mass of the neutralino is much 
smaller than the mass of the stau and Hr > 4, thus VIB is not expected to be relevant. 

It is illustrative to compare these results with other strategies for indirect DM searches. 
In Eigure [3l we show the total neutralino annihilation cross section in the Galactic halo as a 
function of the neutralino mass for the points that fulfil all the experimental limits, including 
the Fermi-LAT constraint from continuum gamma-ray emission from dSphs, currently the 
most stringent bound on {av)o for the range of DM masses considered in this work. As a 
reference, we have depicted the Fermi-LAT upper bounds from dSphs, assuming annihilation 
into T^T~ or hh (solid and dotted lines respectively). Black dots in this plot correspond 
to the few points that exceed the Fermi-LAT bound in Figure [21 For concreteness we have 
chosen the Einasto profile but the results are virtually insensitive to the halo choice. As 
pointed out in Refs. [3I1E51E3H75I, the points with a larger (cju)o generally correspond to 
those with resonant annihilation through the lightest pseudoscalar. Note that these regions 
are significantly less populated in our scan, the main reason being the new bounds on the 
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Figure 2: Thermally averaged neutralino annihilation cross section into two photons in the Galactic 
halo as a function of the neutralino mass. All the points fulfil the experimental constraints, including 
bounds from direct detection experiments and Fermi-LAT data on dSphs, and have a relic abundance 
0.001 < < 0.13. Each colour represents different dominant annihilation channels in the Galactic 

halo, and the different plots correspond to different assumptions on the halo model. 

Higgs sector (very light pseudoscalars can contribute significantly to the invisible SM Higgs 
branching ratio [64(176] ) and the improvement in the determination of low-energy observables 
(with the measurement of BR(H 5 —>■ playing an important role). We should emphasize 

at this point that the resonant condition for annihilation in the DM halo is satished for a 
much smaller mass range than in the early Universe, given the much smaller WIMP velocity 
dispersion in the halo. For this reason, the predicted ^^(cru)o is in general very small. 

Notice that, in principle, one can also have box-shaped features [10311104] in the gamma- 
ray spectrum for hnal states that contain scalar or pseudoscalar Higgses (which can decay 
in-flight to a pair of photons). However, comparing Figures [2] and [3] we observe that the 
total annihilation cross section for is very small (of the order of (cru)o < 10“^® cm^ s“^ 

for points outside the Z resonance), suggesting that the resulting contribution to a box¬ 
shaped feature is several orders of magnitude below the current Fermi-LAT sensitivity for 
line searches. It is also worth noting that variations in the soft slepton mass, m^.g, and 
lepton trilinear coupling, Ae, might lead to points in which the stau mass is similar to that 
of the neutralino. Given the current lower bounds on the stau mass, this is only possible 
for neutralino masses above ~ 100 GeV. Even though we do not expect an increase in the 
contribution to gamma-ray lines, these points might lead to observable features due to VIB 
m- Similarly, a VIB enhancement can occur in the W~^W hnal state [78| when neutralinos 
and charginos are almost degenerate, although it has been argued that this enhancement is 
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Figure 3: Velocity-averaged annihilation cross section of neutralinos in the Galactic halo as a function 
of the neutralino mass. All the experimental constraints, including bounds from direct detection 
experiments and Fermi-LAT dSph data have been taken into account. The solid and dotted lines 
correspond to the Fermi-LAT upper bounds from dSphs, assuming annihilation into or 

respectively. Black dots correspond to those that exceed Fermi-LAT constraints on gamma-ray lines 
in Figure [2] for the Einasto profile. 

more prominent in the limit of large neutralino masses m- 

Let us now put these results in a wider context and contrast them with direct detection 
techniques. In Figure 01 we display the spin-independent neutralino nucleon scattering cross 
section, as a function of its mass. Direct detection constraints are implemented 

using the procedure of Ref. IHZl from SuperCDMS and LUX results. The upper bounds of 
all the direct detection experiments for a Standard Halo Model (SHM) profile and only SI 
interactions are shown for comparison by means of solid lines. As a reference, we also show 
the limits for CRESST [95], and SuperKamiokande [105] (the latter is extracted assuming 
annihilation into Interestingly, the region of low mass neutralinos (with masses as 

small as 3 GeV) is still viable and the predicted for a large fraction of the points 

found lies within the reach of second generation experiments, such as Suj^rCDMS [106] and 
LZ |107j . whose projected sensitivities are shown by means of dashed linefl Notice, however, 
that a substantial region of the explored parameter space lies below the predicted bound 
where coherent neutrino scattering becomes a background for these searches. Our results for 
versus m^o are comparable to those obtained in Ref. |U9] (which were obtained to fit 
the GCE), but span a larger range of solutions. 

^We remind the reader that all the constraints and prospects have been rescaled according to a local DM 
density of 0.4 GeV/cm®. 
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Figure 4: Spin-independent neutralino-proton cross section as a function of the neutralino mass. 
All the experimental constraints, including bounds from direct detection experiments and Fermi-LAT 
data on dSphs have been considered. Solid lines represent the current experimental upper bounds from 
direct detection experiments, whereas dotted lines are the projected sensitivities of next-generation 
detectors. Both of them correspond to a rescaled SHM. Black circles correspond to points whose 
exceeds the Fermi-LAT bounds for the Einasto profile. 

Regarding the results for the spin-dependent contribution with protons, and neu¬ 
trons in general they are rather small, below the current upper bounds derived by 

PICO-2L |108j (which dominates for SD-proton) and LUX (which supersedes XenonlOO |109j 
for SD neutron). We show these results in the upper row of Figure [5l In these plots, we 
also quote published bounds on from the neutrino telescopes, SuperKamiokande [105j 
IceCube mni, BAKSAN m. and ANTARES which produce competitive results for 

DM particles which annihilate preferentially into or W~^W~. Remarkably, a small (but 

non-negligible) part of the scanned parameter space has been excluded based on the sensitiv¬ 
ity of LUX spin-dependent cross section with neutrons. The plot on the bottom left of Figure 
[5] shows the spin-dependent coupling of the lightest neutralino to neutrons, a„, and protons. 
Up. As we can observe, there is a perfect correlation between both quantities (o^ ~ —O.Sop). 
This is a result of the different quark composition of neutrons and protons and the fact that 
in this case the dominant contribution to spin-dependent interactions is the exchange of a Z 
boson (squarks are heavy in our scan). 

Finally, the plot on the bottom right of Figure [3] represents the spin-dependent versus 
spin-independent cross sections (with protons in both cases). This figure can be understood 
as a reanalysis of low-mass neutralino dark matter in the NMSSM with updated constraints. 
Since the direct detection constraints are a function of the DM mass, it is difficult to visualize 
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Figure 5: Top row: spin-dependent neutralino-proton (left) and spin-dependent neutralino-neutron 
(right) cross sections as a function of the neutralino mass. Lower row: a„ and Op coefficients for 
the spin-dependent cross section (left) and spin-dependent versus spin-independent neutralino-proton 
scattering cross section. Solid lines represent the limits from direct detection experiments for a rescaled 
SHM. All the experimental constraints, including bounds from direct detection experiments and Fermi- 
LAT data on dSphs have been considered. No colour code is use to distinguish different annihilation 
channels and black dots represent points that exceed Fermi-LAT constraints on gamma-ray lines in 
Figure [2] for the Einasto profile. 

an exclusion line in this plot. These results perfectly exemplify the necessity of exploring this 
parameter space with complementary targets which are sensitive to both contributions. 

To summarise, the constraints on monochromatic gamma-ray emission for the case of 
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neutralino DM in the NMSSM are in general less important than the dSph bounds and direct 
detection limits. Only in the narrow resonances (mainly with the Z boson) we have found a 
reduced number of points for which exceeds the current constraints. 

3 Spectral features from RH sneutrino annihilation 

The NMSSM can be extended with RH neutrinos in an attempt to accommodate a see-saw 
mechanism that explains the smallness of neutrino masses [531154] . To this aim, in Ref. m 
a gauge singlet RH neutrino superfield, N, that couples to the superfield S was included 
in the NMSSM superpotential. In this construction a Majorana mass term for the RH 
neutrinos is generated after radiative Electroweak symmetry breaking, which is naturally of 
the order of the EW scale. A see-saw mechanism with an EW scale Majorana mass generally 
results in small values for the neutrino Yukawa coupling. This, in turn, implies that in this 
construction the left and right-handed fields have a negligible mixing. It was previously 
shown that the lightest RH sneutrino, Ni, can be a viable DM candidate in wide areas of the 
parameter space |55llll3j , including the low-mass range |9Illll4j , and might also account for 
the GCE mg. 

This model includes five new parameters. Three of them, namely the soft mass parameter 
of the RH sneutrino field, m^y, the coupling between the RH sneutrino and the singlet Higgs, 
Atv, and the corresponding trilinear term in the Lagrangian, A\^, are related to the mass 
and couplings of the RH sneutrino and determine its phenomenological properties as a DM 
candidate. The other two parameters: the neutrino Yukawa coupling and trilinear parameter 
are in general small and have a negligible impact on the DM sector. The RH sneutrino can 
annihilate into a wide range of final states. As in the case of the neutralino, the presence of 
light scalar and pseudoscalar Higgs bosons is very beneficial to reproduce the correct relic 
abundance while fulfilling the LHC constraints on the SM-like Higgs couplings. 

The RH sneutrino in the NMSSM can give rise to a complex spectrum, displaying lines 
and box-shaped spectral featurej§. Regarding the monochromatic gamma-ray emission from 
RH sneutrino annihilation, the only final states that contain photons at the one-loop level are 
77 and Z'j (as the final state violates spin conservation). The RH sneutrino interacts with 
SM particles through the Higgs sector, with a coupling strength controlled by the parameter 
Atv- This leaves two possibilities for diphoton final states in RH sneutrino annihilation: s- 
channel CP-even Higgs, Hf, exchange (with a loop of electrically charged fermions, sfermions 
or bosons), and direct annihilation through a loop of charged Higgses, H^. The Eeynman 
diagrams for these processes are depicted in Figure [6l and for reference, the corresponding 
matrix elements are given in Appendix [Al 

®For a more detailed discussion on these features see Refs. |103II104] . 
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Figure 6: Contributions to the RH sneutrino annihilation into a pair of photons, through a loop 
or mediated by a and a loop of SM and NMSSM particles. 

Concerning box-shaped contributions to the spectrum, in Ref. |115j . we showed how these 
features were obtained from RH sneutrino annihilation into pairs of boosted light scalar and 
pseudoscalar Higgs bosons, which subsequently decay in flight. In order to incorporate bounds 
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Parameter 

Range 

tan /3 

[4,10], [10,20] 

A 

[ 0 . 1 , 0 . 6 ] 

K 

[ 0 . 01 , 0 . 1 ] 

Ax 

[500,1100] 

A^ 

[-50,50] 


[110,250] 

Xn 

[0.07,0.4] 


[-1100,-500] 


[1,150] 


Table 2: Input parameters for the series of scans used in this work in the NMSSM with RH sneutrinos. 
Masses and trilinear parameters are given in GeV. All parameters are defined at the EW scale. 


on these features, we have derived the corresponding limits on the annihilation cross section 
from the Fermi-LAT gamma-ray line bounds. Afterwards we have applied these limits to our 
prediction weighted by the fractional DM density squared, along the box width. 

3.1 Details of the scan and experimental constraints 

We have implemented the RH sneutrino construction in CalcHEP 3.4.3 |116j model files, 
including the one-loop matrix elements for photon final states in RH sneutrino annihilation. 
We define the input parameters at the EW scale and the scan procedure (using micrOMEGAs 
3.6.9 and MultiNest 3.0) and experimental constraints are imMemented in the same way 
as we have described for the case of the neutralino in Section [ 2 dlPI . Due to the smallness of 
the neutrino Yukawa coupling there is no extra contribution from the RH sneutrino sector to 
any of the low-energy observables. 

We have scanned the same ranges of the parameter space as in Refs. [MllllSj . with values 
of the RH sneutrino mass up to 150 GeV. For our choice of soft terms, the lightest neutralino 
becomes the lightest supersymmetric particle for larger values of the RH sneutrino mass. 
Table [5] summarises the values taken for the different input parameters. Gaugino soft masses 
are taken to be Mi = 350 GeV, M 2 = 700 GeV and M 3 = 2100 GeV, thus satisfying the 
Grand Unihcation relation. Slepton and squark soft masses are equal for the three families, 
= 300 GeV, and uIq = rriuc = = 1500 GeV. Trilinear soft terms are chosen 

to be At = 3700 GeV, Af, = 2000 GeV, A^- = —1000 GeV. All these parameters are defined 
at the EW scale. 

®As in the previous section, it should be noted that we do not attempt to fit the GCE. 
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3.2 Results 


In Ref. |91j . we studied the phenomenology of this construction and performed a similar 
scan (limited to lower RH sneutrino masses). From that analysis, we concluded that since 
the new parameters in the RH sneutrino sector can be varied without affecting the Higgs 
masses and couplings, it is much easier to reproduce the correct DM relic density than 
in the neutralino case presented in the previous section. Viable points of the parameter 
space include resonances with the lightest scalar Higgs (we assume no CP violation in the 
sneutrino sector and thus no direct coupling to the pseudoscalar Higgs), but are not limited 
to them. In particular, annihilation into a pair of very light scalar or pseudoscalar Higgses 
is dominant when these channels are kinematically open. We refer the reader to Figures 1 
and 2 of Ref. |91] . where the masses of the lightest scalar and pseudoscalar Higgs bosons are 
represented as a function of the RH sneutrino mass. 

In Figure [71 we show the thermally averaged annihilation cross section of RH sneutrinos 
in the Galactic halo for the di-photon final state, as a function of the RH sneutrino 

mass, for the solutions that satisfy the different experimental constraints mentioned in Sec- 
tion l2.ll Black lines represent the Fermi-LAT bounds at 68% and 95%CL (black dotted and 
thin black lines respectively). Each panel corresponds to a different assumption of the DM 
density prohle, the top left panel corresponds to a NFW, top right to an Einasto, bottom 
left panel to a NEW with adiabatic contraction (NEWc), and bottom right to an isothermal 
density prohle. Since the J-factor for each halo is different the bounds also vary signihcantly. 
Different colours in Eig.jTl represent different dominant RH sneutrino annihilation hnal states 
in the Galactic halo, although we stress that the gamma-ray hux has been computed for each 
point using the contributions for all channels. Thus, bb channel is represented in grey, cc in 
green, in blue, HfH^ in dark blue, in cyan, gg in violet, t+t“ in red, W~^W~ 

in yellow and NN (RH neutrinos) in dark green. Eor clarity, the predictions for each of the 
channels using an Einasto prohle are shown separately in EieurefTOlof AppendixjBj 

As we can observe in Eigure O there is a substantial number of points which present 
an enhanced annihilation into monochromatic photons with respect to the typical value of 
10“^^ cm^s“^, some of them even exceeding the Eermi-LAT bounds. This is due to the 
enhancement of the contribution from the loop diagrams of Eigure [6] under several circum¬ 
stances. On the one hand, these terms are intensihed near a resonance of the RH sneutrino 
with a CP-even Higgs due to the Breit-Wigner effect. This was used in Ref. [l2] in order 
to account for a gamma-ray line at 130 GeV from RH sneutrino annihilation, and in our 
more general scan these resonances occur throughout the whole mass range. In our scan, 
resonances occur with the lightest and second lightest CP-even Higgs bosons. On top of this, 
we have also observed a threshold enhancement [l3] when the particles involved in the loops 
become on-shell. This takes place mainly for loops with W bosons (when mw) and 

charginos (when ^ m-± ~ /u). 
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Figure 7: Thermally averaged RH sneutrino annihilation cross section into two photons in the Galactic 
halo versus the sneutrino mass. All the points fulfil the experimental constraints, including bounds 
from direct detection experiments and Fermi-LAT data on dSphs. Each panel corresponds to the 
Fermi-LAT upper bounds derived for a different DM density profile in an optimised ROI. Colours 
represent different annihilation channels in the Galactic halo, bh final states are displayed in grey, cc 
in green, in blue, in dark blue, A°A° in cyan, gg in yellow, T+r“ in dark red, W~^W~ 

in violet and in dark green. 

All these effects are more prominent in the RH sneutrino than in the neutralino case, 
mainly due to the larger flexibility in the Higgs sector, which is mainly unaffected by variations 
in the mass and couplings of the RH sneutrino (other than through the invisible decay width). 
As a result the predictions for can be larger and a fraction of points in the mass 

range between 30 GeV and 120 GeV, approximately, exceed the Fermi-LAT bound from 
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gamma ray lines search. 

Moreover, the RH sneutrino annihilation into either a pair of very light pseudoscalar or a 
pair of scalar Higgs bosons are dominant annihilation channels when they are kinematically 
allowed [mini]. These Higgses can subsequently decay in-flight into a pair of photons, as 
already mentioned in Ref. |115j . In our scan the annihilation cross section into and 

can be rather large, thereby entailing a sizable contribution to box-shaped features. 
As we explained in Ref. [115| , in order to determine whether these features are observable in 
Fermi-LAT data, we approximate their contribution as a continuum of lines extending from 
the minimum to the maximum box energy. For this reason, excluded points do not necessarily 
correspond to points with a C'^icrv)^^ above the Fermi-LAT bound. In fact, we have observed 
that the number of examples excluded by box-shaped features is larger than those excluded by 
This constitutes another important difference with the neutralino case, for which 
these contributions were more difficult to enhance. 

Finally, we have also found points with a simultaneous enhancement of the Hfgg and 
iL ?77 effective couplings when a singlet-like scalar Higgs with |5'^p/5'j;^p| > 5tan/3 is present. 
In principle, these points might be tested in LHC searches for new light Higgs bosons decaying 
into two photons, since Higgs production through gluon fusion is enhanced with respect to 
the SM expectation. 

Let us now compare these results with the Fermi-LAT bounds from dSphs. Figure[8] 
shows the theoretical predictions for the annihilation cross section of RH sneutrinos in the 
Galactic halo as a function of the RH sneutrino mass for the same choices of DM density 
profiles. All the experimental constraints, including the one from Fermi-LAT dSph data have 
been considered. As a reference, the Fermi-LAT upper bounds from dSphs assuming pure 
annihilation into or bb (solid and dotted lines respectively) are also displayed. Black 

dots correspond to the points of the parameter space for which the monochromatic gamma-ray 
emission and/or the contribution from box-shaped features exceed the Fermi-LAT bounds, 
whereas grey dots correspond to the rest of the points without distinguishing annihilation 
channels. As we can observe, contrary to what happened in the neutralino case, the bound on 
can be more constraining that the bound on the continuum gamma-ray emission for 
a significant region of the parameter space. The amount of excluded points depend slightly 
on the choice for DM halo (NFWc being obviously more constraining). 

Lastly, in Figure [9] we have plotted the theoretical predictions for the spin independent 
RH sneutrino-proton cross section and compared it with current and future direct detection 
bounds (the RH sneutrino, being a scalar field, has no contribution from spin-dependent 
interactions). Once more, the points ruled out by the Fermi-LAT search for spectral features 
are shown in black and each panel represent a different halo profile like in the previous figure. 
It should be stressed that since the RH sneutrino scattering off quarks is mediated by a 
Higgs boson in the t-channel at tree level, in general we expect the RH sneutrino couplings to 
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Figure 8: Thermally averaged RH sneutrino annihilation cross section in the Galactic halo as a 
function of the RH sneutrino mass. All the experimental constraints, including bounds from direct 
detection experiments and Fermi-LAT dSph data have been considered. The solid and dashed lines 
correspond to the upper bounds on (crvjo derived from an analysis of dSph galaxies for pure 
and bb channels, respectively. Black circles correspond to points whose for lines and/or 

box-shaped features exceeds the Fermi-LAT bounds for each DM density profile considered. 

protons, fp, and neutrons, fn, to be very similar, resulting in Notice however 

that the different composition of protons and neutrons in valence quarks, u and d, favours 
the occurrence of accidental cancellations in either coupling when the matrix elements of 
the scalar Higgs boson have different signs (an effect that has been already described for 
the neutralino |117j l. Due to the large flexibility of the Higgs sector in this construction, 
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Figure 9: Spin independent sneutrino-proton cross section as a function of the sneutrino mass. 
All the experimental constraints, including bounds from direct detection experiments and Fermi- 
LAT data on dSphs have been considered. Solid lines represent the current experimental upper 
bounds from direct detection experiments, whereas dotted lines are the projected sensitivities of next- 
generation detectors. Both of them correspond to a rescaled SHM. The dashed line corresponds to an 
approximate band where neutrino coherent scattering with nuclei will begin to limit the sensitivity 
of direct detection experiments. Closed contours represent the areas compatible with the observed 
excesses in DAMA/LIBRA (orange), CRESST (red), CDMS II (blue), and CoGeNT (green). Black 
circles correspond to points whose ^^{crv)^^ for lines and/or box-shaped features exceeds the Fermi- 
LAT bounds for each DM density profile considered. 

such “blind spots” are frequent in our analysis. The ratio fn/fp can be sizable (or very 
small) in such points |1I8| (although the error induced by nuclear uncertainties [I19j is also 
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large). For consistency, as explained above, we are placing direct detection constraints on 
the total nnmber of events, thereby including the contributions from protons and neutrons 
in a consistent way. As we did in the case of the neutralino, the upper bounds for each 
halo have been computed using the procedure detailed in Ref. [971. Remarkably, Fermi-LAT 
bounds on gamma-ray spectral features can rule out many solutions which escape current 
direct detection bounds and are beyond the reach of future experiments. 

To summarise, contrary to what we have observed in the neutralino case, the contribution 
of RH sneutrino annihilation to spectral features in the gamma ray spectrum can exceed the 
Fermi-LAT constraints for a wide region of the parameter space. On the one hand, the Breit- 
Wigner effect near the resonances with CP even Higgses can lead to a sizable annihilation 
cross section into two photons. More importantly, the occurrence of CP even and CP odd 
Higgs final states, and their subsequent decay in-flight gives rise to observable box-shaped 
features. Consequently, the bounds on spectral features can be more constraining than those 
on dSphs and direct DM detection. 


4 Conclusions 

In this article, we have computed the annihilation cross section in the Galactic halo of two 
supersymmetric DM candidates: the neutralino in the NMSSM, and the RH sneutrino in an 
extended version of the NMSSM. For both models, we have studied the resulting gamma-ray 
spectrum, paying especial attention to the contributions from lines (monochromatic gamma- 
rays) and box-shaped spectral features. We have compared the results with the upper bounds 
placed by the Fermi-LAT satellite, using the recent Pass 8 data set. These results have been 
contrasted with the ordinary search for a continuous emission of gamma-rays from dSphs, as 
well as with direct detection experiments. 

To this aim, a series of scans over the parameter space of both models were conducted in 
the region with DM masses below 200 GeV. All the recent experimental constraints have been 
applied. We have implemented the LHC bounds on the masses of supersymmetric particles, 
Higgs mass and couplings, and low-energy observables such as the branching ratios of rare 
processes. In addition, we have incorporated direct detection constraints on the elastic DM- 
nucleus cross section from the SuperCDMS and LUX experiments, computed according to the 
different DM density profiles of the Fermi-LAT ROIs optimised for gamma-ray line searches. 
Finally, we have considered the Fermi-LAT bounds on dSphs. 

In the case of the lightest neutralino in the NMSSM, the majority of the viable points 
of the scanned parameter space occur in the vicinity of the resonances with either a light 
pseudoscalar, a scalar Higgs boson or a Z boson. In agreement with previous analyses, we 
have found that the resulting annihilation cross section into two photons is generally very 
small, beyond the sensitivity of the Fermi-LAT. Only a reduced number of points would 
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be observable thanks to a Breit-Wigner enhancement, which is more effective close to the 
resonance with the Z boson, when m^o ~ 45 GeV, due to its larger decay width. Similarly, 
the contribution from box-shaped features in the spectrum is beyond the reach of Fermi-LAT. 
Even though the predictions of indirect detection through the contribution to the continuum 
gamma-ray flux are also below the current bounds from dwarf spheroidal galaxies, we have 
found that many of these points can be within the sensitivities of future direct detection 
experiments. Additionally, we have provided theoretical predictions for the spin-independent 
and spin-dependent scattering cross section and compared them with the SuperCDMS and 
LUX results. 

On the other hand, in the case of the RH sneutrino in an extended NMSSM, we have found 
that the search for spectral features can be much more constraining. Due to the flexibility 
in the Higgs sector, the RH sneutrino resonant annihilation through a scalar Higgs is much 
easier to achieve than in the case of the neutralino, hence the Breit-Wigner effect near the 
resonance can increase the annihilation cross section into a diphoton final state for a wide 
range of masses. On top of this, RH sneutrino annihilation into light pseudoscalar and scalar 
Higgses, and the subsequent decay in flight of these Higgses into pairs of photons give rise 
to box-shaped features in the spectrum that can be within the reach of Fermi-LAT. As a 
consequence, the search for features in the gamma-ray spectrum by Fermi-LAT can be more 
constraining than Fermi-LAT bounds from dSphs for sneutrino masses in the approximate 
range between 30 and 120 GeV. 

Many of the points that can be probed through the search for spectral features might es¬ 
cape direct detection, thus our findings reinforce the importance of a complementary strategy 
for DM searches. Moreover, the very different predictions for neutralinos and RH sneutri- 
nos demonstrate that the observation of spectral features in the gamma-ray spectrum might 
not just be a smoking gun for DM annihilation, but also could help to discriminate among 
different WIMP candidates. 
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Appendix 


A Matrix elements for RH sneutrino annihilation into two 
photons 


We include here the matrix elements corresponding to the RH sneutrino annihilation into two 
photons. The one-loop diagrams have been shown in Figure [U and we make use of the RH 
sneutrino vertices that can be found in Ref. m- Hereafter, we will express our calculations 
in terms of the function /(r) dehned in Ref. |120j . 


Kt) = 


arcsin T 




log 


1 2 


for r > 1 
for r < 1 


and for a given particle X, we define tx = 
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In the following, we consider rx > 1- 
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where Qt is the top charge and Nc is the flavour number. 
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- Sfermion loop. The leading contributions are expected to come from stops, sbot- 
toms and staus. 
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where the coupling can be found in Ref. |121| . 


- Chargino loop. The squared matrix element can be written as 
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-A ^Hfxtxa "^Xa ^HOxtXb ^xt 
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where ^C'r/o^+y- = —KQiab^L + QibaPR.)- The definition of the Q matrix can be found 
in Ref. |121] . 


B Supplementary plots for RH sneutrino DM 

For clarity, we include in this appendix some complementary plots that help to understand 
how the different annihilation channels are distributed in the parameter space. In particular. 
Figure [To] displays the thermally averaged RH sneutrino annihilation cross section in the DM 
halo, as in Figure [T] but where the points with different dominant annihilation channels are 
presented separately (thus avoiding overlapping). For example, this figure shows that bb is 
widely extended for RH sneutrino masses below 60 GeV, and that and are also 

easy to obtain. 

Finally, in Figure [H] we separate our results in three ranges in the RH sneutrino abun¬ 
dance. If we restrict ourselves to values of compatible with the total DM abundance, 

many of the points below ~ 10 GeV disappear. Otherwise there is no appreciable 

change. Notice that underabundant DM scenarios imply larger values for annihilation cross 
section. However, since we consider a suppression factor the resulting gamma-ray flux 
ends up being slightly suppressed, as we can observe comparing the central and right-hand 
panels. 
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Figure 10: Thermally averaged RH sneutrino annihilation cross section into two photons in the 
Galactic halo versus the sneutrino mass. All the points fulfil the experimental constraints, including 
bounds from direct detection experiments and Fermi-LAT data on dSphs, and have a relic abundance 
0.001 < < 0.13. Each colour represents different dominant annihilation channels. 
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Figure 11: Thermally averaged RH sneutrino annihilation cross section into two photons in the 
Galactic halo versus the sneutrino mass. From left to right, the three panels show different ranges in 
the RH sneutrino relic abundance, the leftmost one being in agreement with Planck results. All the 
experimental constraints are included, together with the bounds from direct detection experiments 
and Fermi-LAT data on dSphs. The colour code is as in Figure [71 
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